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Abstract 

The evolution equation for parton multiplicities in quark and gluon 
jets which takes into account the soft gluon interference is solved nu- 
merically using the initial conditions at threshold. If the A;_L-cutoff Qc 
is lowered towards the hadronic scale Qq of a few hundred MeV, the 
jets are fully resolved into hadrons. Both hadron and jet multiplicities 
in e~^e~ annihilation are well described with a common normalization. 
Evidence is presented within this perturbative approach that the cou- 
pling as{k±) rises by an order of magnitude when approaching the low 
energy region. The ratio of hadron multiplicities in gluon and quark jets 
is found smaller than in previous approximate solutions of the evolution 
equation. 
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1 Introduction 



One of the simplest characteristics of the hadronic final state is the parti- 
cle multiplicity. Following the ideas of a soft hadronization mechanism [||,^ 
the observable hadron multiplicity in a quark or gluon jet is expected to fol- 
low closely that of the final QCD partons in its energy dependence. More 
specifically, within the picture of local parton hadron duality (LPHD) the 
primary parton of energy E generates a parton cascade which is evolved down 
to small scales of a few hundred MeV for the transverse momentum cutoff 
Qq. Then, within the so-called modified leading logarithmic approximation 
(MLLA), which takes into account the leading double logarithmic and next to 
leading single logarithmic terms, a good description of the mean charged par- 
ticle multiplicity in e~^e~ annihilation as a function of the primary cms energy 
^/s = Q = 2E is obtained with a formula 

AO:- (Q) = 2K,,X,i-^, ^) + const. (1) 

Here Ag denotes the multiplicity of partons in a single quark jet of energy E 
and A is the QCD scale. From fits in the energy range (5 = 3... 160 GeV one 
typically finds Qo > A ^ 250 MeV and the normalization K^h ^1-2 
Assuming Kch = \Kaii the number of all hadrons (charged and neutral) is then 
about twice as large as the number of partons {Kaii ~ 2). 

Another characteristic of the final state is the mean multiplicity of jets 
■AfjetiQcQ) for a given resolution parameter Q^, or, in conventional normal- 
ization. He = QllQ"^- Early results at low resolution in 0{a1) played an 
important role as a QCD test and in the determination of the running coupling 
as- The calculation of jet multiplicities to all orders in as became feasible us- 
ing the Durham/fcx algorithm which takes as separation parameter for two 
jets the measure 

1 — cos 9,;,- 

which for small relative angles Qij approximates the rescaled transverse mo- 
mentum {k^lY IQ'^ of the jet of lower momentum with respect to that of higher 
momentum. Then the jet multiplicity MjetiVcQ) refers to all jets with sep- 
aration i/ij > Uc. Jet multiplicities in leading and next-to-leading order of 
ln|/c have been calculated P,p|, p!0| , p!I| and, after matching with the full 0{a1) 
results, a good description of jet multiplicities down to Uc ~ 10~^ has been 



y,,=2uAn{ElE^)^^^ (2) 



obtained |]T2|,|T3|. In these calculations the jet multiplicities are derived in 



absolute normalization {Kaii = 1) from 

Q_ Q 

'Qc A 



A/;et^"(Qc,Q) = 2Ar,(^,^). (3) 
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In this paper we investigate whether both observables, jet multiphcities 
and hadron muhiphcities, can be described in a unified way. Indeed, the 
resummation of the perturbative series has been achieved in both cases by 
using the same type of evolution equation. Also, the Durham/ k± algorithm 
of the jet definition (H) coincides for Qc — > Qo with the prescription k± > Qq 
which is applied in the description of hadrons within the LPHD picture. At a 
first sight there is no smooth transition for Qc Qo from to (|l]) because of 
the two additional constants in Furthermore, various jet observables for 
small Uc deviate from the perturbative predictions (see, for example, review 
||14||), and this is often taken as evidence for nonperturbative hadronization 
effects at the resolution scales Qc ~ 1 — 2 GeV. 

We show that these problems - in the case of the mean jet or particle 
multiplicity - disappear if the underlying coupled evolution equations for quark 
and gluon jets are solved with sufficiently high accuracy; here we perform a 
numerical integration to obtain the exact solutions. 



2 Evolution equation for multiplicities 

We start from the evolution equation for the generating functional of the multi- 
parton final state ||T5|,|^ which, in the parton splitting process A B-\-C, takes 
into account angular ordering, energy conservation and the running coupling 
at the one-loop order. By appropriate differentiation we obtain the evolution 
equations for the mean parton multiplicities J\fq and Afg in quark and gluon 
jets with jet virtuality k (see also |Tl|) or with 



r/ = ln^, K = gsin®. (4) 

where denotes the maximum angle between the outgoing partons B and C. 
At fixed cutoff Qc these equations read 

dAfM _ /^-^^^^^£M[$-.(,){Ar,(^ + ln^)+Ar,(r/ + ln(l-;.))-Ar,(r/)} 



dr] Jzc 2n ^ 

+nf^gq{z){Xq{r] + \n z) + Afq{r] + \n{l ~ z)) - ^g{r])}] (5) 



dz^^<l>qg{z){Afg{r^ + \nz)+ Nq{r, + ln(l - z)) - MM)- 
df] Jzc 27r 

The splitting functions ^ab for parton splitting A ^ B |T^ are taken with 
normalization as in 0, i.e., <^gg{z) ~ ^Nc/z at small z {Nc and Uf denote 
the number of colours and fiavours). Because of its symmetry property it 
is convenient to replace ^^gg{z) by the asymmetric function $^^^(2;) = (1 — 
z)%giz) in (§) [0. 
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The boundaries of the integral over the parton momentum fractions z are 
determined by the lower cutoff in the transverse momentum measure defined 
according to (|^) with the approximation zj^"^ ^ — zb 

k± = min(z, 1 — z)k > Q^. (6) 

The lower bound z^ in is obtained for minimal z = Qc/{Q sin ^) which 
is found for G ~ | (in this configuration, because of transverse momentum 
conservation, the high energy secondary parton has production angle "i?^ ~ 
and the soft particle {}b < f); niore generally, 6 corresponds to the half 
opening angle of the jet {-ds < ©)• Therefore Zc is given by 

Zc = = y2yc = e ^ (7) 

Since Zc < |, one finds 1/c < | and r] >\n2. 

The coupling is given by as{k±,nf) = 2Tr/ {bln{k±/ A)) with b = {llNc — 
2n/)/3. We evolve as(A;_L, n/)~^ with k± by matching smoothly the couplings 
for n f and nj + l at twice the heavy quark thresholds according to the formula 



Q{kt - 2mn.) 



Uf 



The factor two in the threshold for heavy quark production takes into account 
that heavy quarks are produced in pairs by a gluon. 

The differential equations are defined only for r] > In 2. The initial condi- 
tions for the solutions of (|^) then read (for any fixed Qc) 

Ug{r]) = Uq{r]) = 1 for < < ln2. (9) 

Analytical solutions using the boundary conditions at threshold have been 
derived within the MLLA for gluon jetsfl^ and for the coupled system of quark 



and gluon jets|TO[. Whereas the MLLA yields a good high energy behaviour, 
it leads to inconsistencies [M'^irj) < 0) near threshold. In the present analysis 
we solve the equations (^ and (j^ numerically^ 

The two jets in e~^e~ annihilation evolve independently in this approxima- 
tion. Near threshold this factorization of the generating functional does not 

^We start from J\fa{rj = In 2) = 1, calculate Afaiv + ^rf) in steps of length 5rj from 
the derivatives (^ as J^^{r])Sri using the trapezoidal rule for the integration and linear 
interpolation for the required M values under the integral. The integration is performed 
with logarithmic variable y = Inz. 
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hold any more and non-logarithmic terms become important. For the k± res- 
olution criterion (|^) the inelastic threshold is found for symmetric 3-jet events 
(9 = 2; = |) with = I or 7] = I In |. This is to be compared with ?/c = | 
and ?7 = ln2 for eq. The difference of thresholds will clearly affect the 
results in the large Uc region. 

An improvement in this region can be achieved by replacing the contribu- 
tion of 0{as) in (|^) by the explicit result for e'^e^ 3 partons in 0{as) (see 
also Ref. |TO|,p|). The lowest order contribution J\f^^^ of the evolution equation 
is obtained by taking the first iteration of (|^) with the initial conditions @, 
i.e., replacing the three A/'-terms in the curly brackets in (^ by unity. The full 
0{as) contribution is found by numerical integration of 

CFas{k±) zf + z? 



Af'-^^'iVc) =2 f dz, r dZ2e{d23 - Vc) 
Jl/2 Jl-zi 
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271 {1-Zi){l-Z2) ^^Q^ 



^23 = min ( — , — ) (1 - > ?/c (11) 

\Z3 Z2/ 

where Zi (2:2) denote the quark (antiquark) and z^ = 2 — zi — Z2 the gluon 
momentum fractions [Cp = 4/3). Here the coupling is taken again as running 
with the k_i_ measure (0) kj_ = (c?23(5^)2 as scale in the integration region 
Z2 < z\ in (10). The 0[a^ corrected solution M^^^^ is then obtained from 

AC:r(yc) = 2Ar,(y,) - 2M^^\y,) +U'-^^\y;). (12) 

A corresponding improved result for the gluon jet is not yet available and 
has to be considered for each process separately. In this paper we apply the 
equation (|^) to the full range of scale parameters Qo < Qc < Q, i-e., we study 
the jet region (y^ > 0.01) and also the transition from the jets to the fully 
resolved hadrons (?/c ~ {Qo/QY)- 



3 Confronting predictions with experiment 
3.1 Multiplicities in e+e annihilation 

In Fig. la we show the data of the average jet multiplicity at Q = 91 GeV 



T2| , p^ as a function of the resolution parameter Uc = Qc/Q^ as obtained 
with the k± algorithm. The theoretical predictions from (|1^) for the jet data 
are given in absolute normalization in terms of the single parameter A. Also 
shown are the data on hadron multiplicities in the energy range Q = 1.6. . . 
91 GeV taken as Afaii = ^J^ch as a function of the same scale parameter, 



now calculated as yc = Ql/Q"^, where the parameter Qq corresponds to the 
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parton k± cutoff characterizing a hadronic scale and is obtained from a fit to 
tfie data. Another adjustable parameter here is the overall normalization Kaii 
which relates the parton and hadron multiplicities in Afaii = KaiiAfcorr ■ We 
determine first the A parameter from the jet multiplicity (lower data set in 
Fig. 1) and then the Kau and Qq from the hadronic multiplicity (upper data 
set). A good description of the data is obtained with parameters 

K^ii = 1, A = 500 ± 50 MeV A = In ^ = 0.015 ± 0.005 

^ (13) 

which correspond to the curves shown in Fig. la. The quantity M — 2 derived 
from the same experimental data is shown in Fig. lb. The dashed curves 
represent the contribution 2Mq in ([T2|) from the evolution equation alone. The 
fits can be seen to describe quantitatively the hadron multiplicity, the jet 
multiplicity for the small and large ?/c, whereas in the intermediate region 
around ijc — 0.01 the fit exceeds the data on A/" — 2 by up to 20% (up to 5% in 
M). Deviations of this type may indicate the importance of two-loop 0{a1) 
terms not included in the analysis. 

As a remarkable result of our analysis we find that the common normaliza- 
tion Kail = 1 is possible without difficulty; the normalization parameter Kaii 
is correlated with Qq and can be varied within about 30%. The solution with 
Kail = 1 is natural as it provides the correct boundary conditions A/" = 2 at 
threshold for both hadrons and jets. Then a continuous connection between jet 
and hadron data results, since both are described by the same evolution equa- 
tion (|12D. In this description a single hadron corresponds to a single parton of 
low virtuality Qq. 

It is interesting to compare the two sets of data at the same Uc- The differ- 
ence between the two curves comes entirely from the running of as- Namely, 
if as were kept fixed in equations (|^) the resulting multiplicities would depend 
only on the ratio of the available scales through the variable rj = — |ln(2?/c) 
but not on the absolute size of Q and the two curves in Fig. 1 would coincide. 
A model with fixed ag would predict at high energy a power like dependence 
on Q, i.e. a straight line in Fig. la between the two curves shown. 

The largest difference occurs near the inelastic threshold at j/c ~ | or 
r/ ~ ln2. In this region A/" — 2 is dominated by the contribution of order 
and the jet and hadron results at the same i/c are in the ratio of the relevant 
typical coupling constants. At the lowest available energy for the hadron mul- 
tiplicity (?/c ~ 0.1, Q = 1.6 GeV) this ratio of couplings is larger than 10, 
so the typical coupling to produce hadrons at 1.6 GeV is as > 1. The good 
matching of the prediction for hadrons with the boundary condition at thresh- 
old suggests that the coupling is rising even more towards lower energies. The 
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strong variation of the coupling at small scales has been found important also 
in the description of particle energy spectra|]5|,p^. On the other hand, the 
jet multiplicity is first rising very slowly with decreasing because of the 
small coupling (a^ ~ 0.1). Only for very high resolution, if is lowered from 
900 MeV ~ 10""^) to the final 500 MeV, about three quarters of the final 
multiplicity are produced, three times as many jets as in the large complemen- 
tary kinematic region. The steep rise toward small Uc in Fig. 1 reflects the 
singularity in the coupling at |/a = A^/Q^, which is however screened by the 
hadronization scale Uq = {QoY/Q'^ > Va- This behaviour is qualitatively de- 
scribed by the high energy DLA resuh[|] J\f ~ Ko{Ay/Tc) h{^\/M^^J^) + ■■■ 



with A = \JlQNc/h which at high energies and small Ac = ln((5c/A) behaves 

as In Acexp(A-y/ln(K/A)). This form, for fixed Qc, describes the slow rise of 
multiplicity with k, ^ Q, whereas, for fixed Q, it exhibits the logarithmic sin- 
gularity for Ac — * 0. This singularity arises in the first iteration of (the 
"Born-term") from the contribution Jq dy/{y + A). 

Previous approximate solutions had been restricted to i/c > 0.002 for jets, 
whereas for hadrons the analytical results could only be applied to the higher 
energy region with a larger Kch factor.^ Our results are consistent with the 
previous finding A < 0.1 from the energy moments on the other hand, A = 
is not allowed in (^), as the multiplicity would diverge in this case. 

It will be interesting to study this behaviour further near the transition 
Qc Qo also at other cms energies, especially for the jet multiplicity, since it 
varies strongly in the low i/c range. In the experimental analysis using eq. @ 
the full hadronic multiplicity will only be reached for i/c — > 0. The same lower 
cutoff for partons and hadrons is achieved if we interpret Qc as transverse 
mass cutoff with Qo as effective mass parameter (see also P,PD|), i.e., if we 
relate the theoretical predictions and experimental results from cutoff (0) by 

yf = yT' + Ql/Q'- 

The value ( p!3D for the scale A is larger than that found in previous analyses 
of particle spectra[P,§J^. This difference comes from taking the Durham k± not 
only in the integral boundaries but also as argument of and in the definition 
of hadrons by k± > Qq. We performed an alternative calculation where we 
used the standard k± < k± as argument of as (i.e., k± = z{l — z)n instead 
of eq. (P) and the exact k± in eq. (10)) and inserted the cutoff Q{k± — Qq) 
into eqs. (|^) and (10); then a description of comparable quality is obtained 

^The different value of the normalization factor as compared to previous analyses (for 
example ^,^,||) is due to the use in the MLLA approximation of the initial conditions which 
have negative slope A/"' (77 = 0) < and therefore lead to a delayed rise with energy pO[ |. 
Our exact solution in the present cms energy range is about twice as large and rises more 
slowly than the analytical MLLA solution for Qq — A p7|. 
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with the lower scale A ~ 0.35 GeV and the same A parametei]^. The threshold 
for hadrons is shifted in this case to the lower value Uc = 1/12 and also the 
hadron data points are shifted to the left in Fig. 1 according to the lower value 
of Qq. Our conclusions about the order of magnitude variation of remain 
unaltered. The ambiguity in choosing the scale for «<, may be reduced in 2- 
loop calculations; for now we stay with the conceptually simplest possibility 
and take k± throughout. 



3.2 Multiplicities in gluon jets 

The evolution equations (|^) yield also results on the mean multiplicity 
in gluon jets. The results derived here refer to the multiplicity in the full 
hemisphere of the gluon jet and can be measured, for example, in a final 
state with a primary gg colour singlet state. Such a state is expected to be 
produced in the decay of heavy quarkonia; recently the multiplicity in the decay 
^ ~^ 991 has been measured by the CLEO Collaboration [^. An initial state 
of this type can also be realized approximately in e^e" 3 jets with nearly 
antiparallel q and q recoiling against the gluon p|,^; this configuration has 



been analysed by the OPAL Collaboration [^0. Our results do not apply 
immediately to other configurations (like "Y" or "Mercedes") where the large 
angle soft radiation from the gluon is suppressed (for reviews, see |2^,|T^). 

Of special interest is the ratio of gluon and quark jet multiplicity r = 
Ng/Mq, which is predicted to approach asymptotically the value r = | 



The next-to- leading order corrections decrease this ratio at finite energies [25 



and yet further with the inclusion of higher order corrections [26,^|. There is 



an uncertainty with these asymptotic predictions in that the region of validity 
is not known a priori. A complete solution of the evolution equations (|^) 
requires an initial value for the multiplicities at some energy. Such results with 
absolute normalization at threshold have been presented within the MLLA . 



We have obtained the results from our evolution equation (^) with initial 
condition (^), however, we have not yet included a low energy correction term 
as in ( |T^ ) which is beyond the scope of this paper. The results for the ratio r 
are shown in Fig. 2 again as function of Uc = {Qc/QY for hadron multiplicities 
in comparison with the OPAL and CLEO data mentioned above. Also shown 
is the prediction for jets at fixed Q = 91 GeV which could be tested using the 
same type of events with antiparallel qq. 



■^The remaining difference to the earlier A ~ 250 GeV comes from taking the scale 

y/2E in (I) instead of k = i?. 
''The full multiplicity in the hemisphere opposite to qq corresponding to our calculation 
is obtained if the higher jet multiplicities are also taken into account. 
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There is a very good agreement of our calculation with the OPAL data 
which refer to jet energies of ~40 GeV. Prom Pig. lb one may expect that 
the inclusion of the full 0{as) corrections would not modify the prediction 
essentially (by more than about 10%). On the other hand, the CLEO data 
fall considerably below the prediction r ~ 1.25 from the evolution equation. 
In this region, however, the low energy corrections may not be small: in Pig 
lb these corrections correspond to a factor of two. A more precise prediction 
of the Y data will require the inclusion of the full 0{as) correction to the 
Y —>■ gg'j process. 

4 Conclusions 

We solved the QCD evolution equations for quark and gluon jets exactly by 
numerical integration and also supplemented the full 0{as) correction for e'^e~ 
annihilation. With this improved accuracy beyond MLLA we obtain a unified 
description with common absolute normalization of the mean jet multiplicity 
at LEP-1 and the hadron multiplicities in the energy range from 1.6 to 91 
GeV within the errors or at least within 5%. In this description only the QCD 
scale A and a k± cutoff parameter Qq > A for the hadrons according to the 
LPHD picture have been adjusted to the data. The ratio of gluon to quark 
multiplicity at LEP-1 is described as well. 

Of particular interest is the success of this model in the regions where 
perturbative QCD is not expected to be relevant a priori, namely for small 
transverse momenta, where the running coupling gets large, either for the 
hadron multiplicity near threshold or for the jet multiplicity at very high res- 
olution. In both cases the strong effect of the running is clearly seen from 
the comparison of the two multiplicities at the same scale yc] taking the trans- 
verse momentum as the argument of as provides about the right separation. It 
appears that perturbative QCD is still successful in these extreme kinematic 
regions for this inclusive quantity. 

The model for hadron production emerging from this analysis is very sim- 
ple. A hadron can be treated formally like a narrow microjet of partons but 
with the radiated partons being so soft (A < A;_l < Qo) that they cannot be 
resolved and are therefore confined in a region characterized by the hadronic 
scale Qo of a few hundred McV. So a hadron in this picture docs not cor- 
respond to a colour singlet state of "valence partons" but rather to a single 
parton (quark or gluon) with about the same momentum together with a "sea" 
of soft confined partons which take care of the colour neutrahzation. In this 
way the soft hadronization picture can become consistent with the perturbative 
treatment of particle production. 
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Figure Captions 



Fig. la: Data on the average jet multiplicity at Q = 91 GeV |[T^ . [13[| and the 



average hadron multiplicity (assuming A/" = ^Mch) at different cms energies[|8 
with Qc = 0.507 GeV as a function of ye- The solid (dashed) line shows the 
prediction in absolute normalization for the hadron (jet) multiplicity obtained 
by using eq. (|12|) with parameters from eq. (p^. The right most data point 
for hadrons {Qq = 1.6 GeV) refers to pions only. 

Fig. lb: Data as in Fig. la, but A/" — 2 is now shown. The solid lines show 
the complete predictions of eq. (|12D for both the particle and jet multiplicities, 
the dashed lines the contribution of eq. (^) alone. The same parameters as in 
Fig. la are used. 

Fig. 2: Data on the ratio of average hadron multiplicities in quark and gluon 
jets at ?/c ~ 3 X 10"^ from 0PAL|2g and around yc ~ 0.01 from CLEO[|g 



with Qc = 0.507 GeV as a function of yc- The solid (dashed) line shows the 
prediction for the ratio of hadron (jet) multiplicities in quark and gluon jets 
obtained by using eq. (^ with the same parameters as in Fig. 1. 
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Figure 1: a 
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Figure 1: b 
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Figure 2: 
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